Supermassive recoil velocities for binary black-hole mergers with antialigned spins 
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Recent calculations of the recoil velocity in binary black hole mergers have found the kick velocity 
to be of the order of a few hundred km/s in the case of non-spinning binaries and about 500km/s 
in the case of spinning configurations, and have lead to predictions of a maximum kick of up to 
1300 km/s. We test these predictions and demonstrate that kick velocities of at least 2500 km/s are 
possible for equal-mass binaries with anti-aligned spins in the orbital plane. Kicks of that magnitude 
are likely to have significant repercussions for models of black-hole formation, the population of 
intergalactic black holes and the structure of host galaxies. 
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Introduction. — A well known phenomenon of general 
relativity is the recoil or rocket effect due to the emission 
of anisotropic gravitational radiation 0, 0, Q . The loss 
of linear momentum radiated away in the form of gravi- 
tational waves imparts a recoil or kick on the remaining 
system which then moves relative to its original centre- 
of-mass frame. This effect is particularly pronounced for 
the inspiral and merger of two compact objects and thus 
may have dramatic consequences for the merger of mas- 
sive black holes residing at the centers of galaxies when 
their hosts undergo merger. 

Massive black holes with masses 10 5 to 10 9 ' 5 M© are 
not only known to exist at the centre of many galaxies, 
but also to have a substantial impact on the structure 
and formation of their host galaxies, as is demonstrated 
by the correlation of the black-hole mass with the bulge 
mass, luminosity and velocity dispersion 0, 0, 0| ■ The 
merger of galaxies is thus likely to be accompanied by an 
inspiral of the central black holes. In order to assess the 
impact of the resulting black-hole recoil on questions such 
as the bulge structure, the formation history of massive 
black holes and interstellar and intergalactic populations 
of black holes, it is vital to have a good understanding 
of the kick magnitude and, in particular, the maximum 
possible kick velocities. 

The first study of the recoil effect for inspiralling bi- 
naries was performed by Fitchett Q in the framework of 
non-spinning point particles subject to Newtonian grav- 
ity. The resulting recoil calculated using the lowest or- 
der multipoles is given in his Eq. (3.21) and predicts 
kick velocities of the order of 1000 km/s that exceed the 
escape velocities from massive galaxies The prob- 
lem was reinvestigated using the particle app roximation 
13, EH, El El , post-Newtonian methodsEJEl El, E3] , 
the Close-Limit Approximation 18, 13, 2fjj] as well as nu- 
merical simulations jU [12] • The picture that emerged 
from these studies is that the recoil from unequal-mass, 
non-spinning binaries is unlikely to exceed a few hundred 
km/s. 

Recent breakthroughs in the numerical simulation of 
black-hole binaries [23, 24, [25| have enabled investiga- 



proximations other than the numerical differencing of 
the Einstein equations. First studies addressed the re- 
coil from non-spinning binaries and confirmed the rel- 
atively small magnitude of the kick velocities for mass 
ratios in the range 1:1 to 1:2 [H, [53]. In a previous 
publication [28| , we have presented the most comprehen- 
sive study of this problem and found a maximum kick of 
175.7±11 km/s for a mass ratio 77 = m\rciil (m\ +m2) 2 = 
0.195 ± 0.005. 

More recent numerical studies have shown, however, 
that as expected significantly larger kicks are realized if 
one allows at least one black hole to spin. Simulations 
of equal-mass binaries with spins orthogonal to the or- 
bital plane predict kick velocities of 475 and 440 km/s, 
resp ectively, in the limit of extreme Kerr black holes 
[23, H3]. Kicks of tens of km/s have been obtained from 
head-on collisions of spinning black holes [3l|. These 
results are consistent with the effective-one-body post- 
Newtonian (PN) study in (32J. In addition, Campanelli, 
et. al. [33] obtain Ukick = 454 km/s in the case of a non- 
spinning black hole orbiting a spinning counterpart of 
twice the mass with the spin oriented at —45° relative to 
the orbital plane. They further estimate that an equal- 
mass binary with spins aligned in the orbital plane could 
produce a maximum kick of around 1300 km/s based on 
a PN calculation by Kidder (l5j . 

In this letter we inv estig ate the scenario suggested by 
Campanelli, et. al. fl5L [33I] and find that a kick of 
2500 km/s is possible — larger than the escape velocity 
of about 2000 km/s of giant elliptical galaxies. 

Numerical framework. — The numerical simulations 
presented in this work were performed independently 
with the Bam and the Lean codes. These codes em- 
ploy the moving-puncture method [2^, [24[ , and are de- 



scribed in detail in 34, 3 



We note 

lowing modifications of the Cactus [3 



tions of the recoil problem without any restrictive ap- 



however, the fol- 
and Carpet 13711 
based Lean code relative to the version presented in [351] : 
time evolution of the TwoPuncture initial data [38| is 
performed using the fourth-order accurate Runge-Kutta 
method, the variable (f> has been replaced by the new 
variable x — e an d the code uses fifth-order prolon- 
gation in space. The two codes represent independent 
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v kick (fcm/s) 


MI 


±3.257 


±0.133 


0.363 


±0.2 


2450 ± 250 


Mil 


±4.0 


±0.1125 


0.287 


±0.2 


2650 



TABLE I: Initial puncture parameters and final kick velocity. 



implementations of similar techniques, thereby allowing 
important cross-validation. 

The gravitational waves emitted by the binary are ex- 
tracted by computing the Newman-Penrose scalar 
over spheres of constant radius at different distances. 
The notation and equations used for this procedure are 
explained in detail in section III. A of [34J . In order to es- 
timate the recoil velocity of the system, we compute the 
total linear momentum radiated during the simulation 
directly from ^4 using equation (1) in (28|. The pulse 
of spurious radiation from the initial data, which would 
lead to an over-estimation of the final kick by about 5%, 
is not included in this calculation. We have neglected 
any recoil that may have accumulated during the ear- 
lier inspiral of the black holes; however, for the final kick 
values we obtain, this introduces an uncertainty of less 
than 1% in our results (as obtained from PN estimates 
in the nonspinning case (l6|). which is far less than the 
uncertainty due to other errors. 

We consider configurations with spins of equal size but 
opposite direction lying in the orbital plane. The initial 
parameters of Model I (MI) were chosen without partic- 
ular regard to obtain a quasi-circular configuration. On 
the other hand, the initial momenta used for Model II 
(Mil) were computed using formula (4.7) from [151 ] to 
generate a quasi-circular orbit. The parameters and the 
values of the kick obtained for the two models are sum- 
marized in Table [U The initial configurations, consisting 
of two equal-mass black holes with spins perpendicular 
to the orbital angular momentum, are chosen to maxi- 
mize Eq. (3.31b) in [l5| as suggested in (33[. Each black 
hole has a total mass of m ~ 0.5 and we consider spin 
parameters of the order of a/m ~ 0.8, which is below the 
recent estimates presented in (39j . 

Results. — We performed evolutions of model MI 
with the Lean code using different resolutions and com- 
puted the radiated linear momentum at different ex- 
traction radii. In Figure [T] we present the recoil speed 
v z = —P z /m of the final black hole as a function of time 
for the resolutions 1/36, 1/44 and 1/48. In this context 
we note that the x and y components of v are smaller 
than 1 km/s, so that the z-component is practically iden- 
tical with the total recoil. The bottom panel of the figure 
demonstrates fourth-order convergence of the recoil, and 
we estimate the uncertainty as 43 km/s or 1.5 %. 

The dependence of the recoil velocity on the extrac- 
tion radius is presented in Figure O We find the re- 



sulting error to be reasonably well approximated by a 
l/ r cx falloff and thus obtain uncertainties due to finite 
extraction radii of 120 km/s or 4.5%. Using a conser- 
vative error estimate, the magnitude of the final kick is 
2450 ±250 km/s. 
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FIG. 1: Upper panel: the z-component of the recoil for model 
MI as a function of time for grid resolutions h\ — 1/36, /12 = 
1/44 and hz = 1/48. Lower panel: the differences scaled for 
fourth-order convergence. 
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FIG. 2: Recoil for model MI as a function of time for different 
extraction radii. Extrapolation to r — > 00 yields a magnitude 
of 2450 km/s with an error estimate of 4.5%. Also shown is 
the kick from model Mil using Bam. 

In order to obtain an independent estimate of the kick 
magnitude for this type of spin-configuration, we evolve 
model Mil using the Bam code. This model starts at 
larger initial separation and represents a quasi-circular 
configuration. These simulations thus enable us to assess 
the uncertainties arising from finite black hole separation 
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and deviations from circularity. 

The trajectories of the black holes are shown for model 
Mil in Figure [3] and demonstrate how the black holes 
move out of the initial orbital plane and acquire substan- 
tial momentum in the z-direction after the merger. 
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FIG. 3: Coordinate positions of the black-hole punctures for 
model Mil up to t = 180. The black holes move out of the 
original plane and after merger the final black hole receives a 
kick in the negative ^-direction. 

In summary, our simulations predict recoil velocities 
of 2500 km/s with a conservative error estimate of 10% 
or 250 km/s. This value can be compared with the 
escape velocities from dwarf elliptical and spheroidal 
galaxies (<■ 300km/s) and from giant elliptical galaxies 
(< 2000km/s) 0. 

Our prediction is larger than any previous numerical 
result by a factor of five, but other numerical studies 
considered quite different initial configurations. Since 
the estimate of 1000 km/s [33| for the spins studied here 
involves a post-Newtonian argument, it is perhaps not 
surprising that a fully relativistic simulation obtains a 
different result. This result is very intriguing, but we 
point out that so far only the two data points presented 
here are available. It will be important to perform studies 
at higher numerical resolution, and to study neighboring 
data sets, for example to vary the initial separation, spin 
configuration and mass ratio of the black holes. 

Astrophysical relevance. — Before we discuss in more 
detail astrophysical implications of our findings, we em- 
phasize one additional important point: while our results 
demonstrate that kicks as large as 2, 500 km/s are possi- 
ble provided the inspiralling black holes have appropriate 
spin alignment and magnitude, it does not enable us to 
make any statements on the probability that these kicks 
are realized in typical astrophysical merger scenarios. A 
conclusive assessment of the astrophysical implications 
therefore requires systematic parameter studies as men- 
tioned in the previous section. 

Recoil velocities close to or exceeding the escape ve- 



locities of the black hole's host give rise to a population 
of black holes away from the galactic nuclei 4l| . In con- 
trast to the non-spinning case, the recoil velocities calcu- 
lated from our simulations exceed the escape velocities of 
even massive galaxies and thus imply larger populations 
of wandering black holes. Observational consequences of 
wandering black holes have been discussed in detail in 
Volonteri and Perna [42] • 

Black-hole recoil has also been found to give rise to a 
core formation in the central density profile in the host 
stellar bulge. Boylan-Kolchin et al. [43j find this effect to 
be most pronounced for recoil velocities just below the 
escape velocity. 

Libeskind, et. al [44j have further found that grav- 
itational recoil manifests itself as a scatter of the rela- 
tion between black hole and bulge mass. They find this 
effect to be sensitive to the stability of the disc under 
galaxy merger. Their simulations indicate a constraint 
of Kkick ^ 500 km/s. Merritt, et. al. [45| find constraints 
of a similar magnitude from narrow emission-line analysis 
of quasar spectra. 

Ejection of black holes affects the rate of binary black 
hole inspirals in globular clusters and thus the predicted 
event rates for gravitational wave observatories LIGO, 
GEO600, TAMA and VIRGO. The latest kick estimates 
for non-spinning binaries have been taken into account 
by O'Leary, et. al. [46|. It will be of interest to estimate 
the impact of spins on such simulations. 

The formation history of massive black holes at z ~ 6 
is often described in the context of hierarchical structure 
formation via accretion and merger of black holes resid- 
ing in dark matter halos (see, e.g. [Ilj]). The ejection of 
black holes from their hosts via gravitational recoil puts 
constraints on the maximum redshift at which the pro- 
genitor seed holes might have started merging and might 
necessitate accretion growth above the Eddington limit 
(ifl 47 1 and may also lead to a population of interme- 



diate mass black holes [42J. Black hole mergers would 



not occur if seed black holes are rare at high redshifts, so 
that black-hole binaries would not commonly form as a 
consequence of dark matter halo mergers [481 ]. 

Finally, black hole recoil has been suggested to mani- 
fest itself directly in observations of off-centre radio-loud 
active galactic nuclei [48j j and off-nuclear ultraluminous 
X-ray sources in nearby galaxies (49l |. 

Astrophysical studies of the consequences of kick ve- 
locities on these scenarios have so far commonly assumed 
recoil velocities significantly below the values obtained in 
this work. It will be interesting to see how these conse- 
quences are affected by our results. It will also be im- 
portant to estimate the kick magnitudes in a far wider 
volume of parameter space. 

Note added after preparation of this manuscript: After 
the first publication on gr-qc of the results reported here, 
there have been two independent computations of simi- 
larly large kicks, first reported in (Hoj and in an updated 



4 



version of 13311 . Since our work is based on the original 
version of [33], it is relevant to note that the authors 
have now added new numerical data to their preprint for 
a configuration showing a kick of 1830 km/s. They have 
also modified their argument for the extrapolation from 
post-Newtonian results (see also [HT[ ) now obtaining an 
upper limit of 4000 km/s [H^]. This is consistent with the 
results we report here and supports the main point we 
make, namely that kicks much larger than previously ex- 
pected are indeed realized in numerical simulations with 
potentially great significance for astrophysics. 
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